June 2011 
SNUTPll-004 



Eigenvalue Distributions in Matrix Models for 
Chern-Simons-matter Theories 



Takao Suyama 







BK-21 Frontier Research Physics Division 
and 

Center for Theoretical Physics, 

Seoul National University, 
Seoul 151-747 Korea 



Abstract 

The eigenvalue distribution is investigated for matrix models related via the localization to Chern- 
Simons-matter theories. An integral representation of the planar resolvent is used to derive the positions 
of the branch points of the planar resolvent in the large 't Hooft coupling limit. Various known exact 
results on eigenvalue distributions and the expectation value of Wilson loops arc reproduced. 



^e-mail address : suyama@phya.snu.ac.kr 



1 



1 Introduction 



Superconformal Chern-Simons-matter theories have attracted much attention recently in relation to the 
dynamics of M2-branes. The worldvolume theory on M2-branes was proposed in [1] [2] [3] [4] [5] based on 
a 3-algebra structure. The theory turned out to be equivalent to a Chern-Simons theory with a product 
gauge group coupled to matters [6]. Chcrn-Simons- matter theories with product gauge groups were also 
studied in [7] in a different context, and then the construction of [7] was extended to more general theories 
[8] (see also [9]). Soon after those developments, the worldvolume theory on any number TV of M2-branes 
placed at the tip of an orbifold was constructed in [T^. Since this theory, known as ABJM theory, allows 
us to take the large N limit, one may discuss its gravity dual, as in the case of A/" = 4 super Yang-Mills 
theory in four dimensions However, as usual in AdS/CFT correspondence, it is very difficult to 

check the correspondence between the field theory and the dual gravity or M-tlieory, unless it is possible 
to perform an exact calculation in the field theory side, such as a calculation of a quantity which is 
protected by the BPS nature or the integrability. 

Fortunately, there exists a technique, known as the localization, which is quite powerful for calculating 
some class of observables in field theories. The localization technique may be used to calculate expectation 
values of operators which preserve a particular supercharge. The reason that this technique is so powerful 
is that it enables one to perform a path-integral exactly by just evaluating the saddle-points. In the context 
of AdS5/CFT4 correspondence, the localization technique was used in [12] to calculate the expectation 
value of the half-BPS Wilson loop. As a result, one obtains a function of the 't Hooft coupling A which 
reproduces, via the relation [E] [14] [15] , the behavior of the string worldsheet in AdSs exactly in the large 
A limit. 

A similar technique was applied to Chcrn-Simons- matter theories in |16j . The localization formula 
obtained in |16| is much simpler than that for a generic four-dimensional Af = 2 gauge theory in |12| . 
Therefore, the analysis using this formula can be done rather easily. First, the ordinary 't Hooft limit of 
ABJM theory was investigated [I7][IH] using the localization formula, and the behavior of the BPS Wilson 
loop [18] constructed in [19][20][2T] and the free energy [22], expected from AdS/CFT correspondence, 
were reproduced exactly. Note that the dual theory for the 't Hooft limit is Type HA theory [lU]. Next, 
another limit, called the M-thcory limit, was investigated in [23]. This is a large N limit with k kept 
finite, so that the dual theory should be M-theory, not Type HA theory. The technique developed in [53] 
enabled the authors to derive the exact behavior of the free energy and the Wilson loops in the M-theory 
limit. Remarkably, this technique can be applied to various other Chern-Simons-matter theories whose 
planar resolvents were not known. 

The investigation based on the localization formula had been limited to the theories which have at least 
A/" = 3 supersymmetry, since otherwise the dimensions of the fundamental fields may be renormalized in 
the IR, and localization formula in [16j may not valid. Recently, the localization formula was generalized 
to A/" = 2 theories in [24] [25] which plays a central role in recent developments on Chern-Simons-matter 
theories, especially on the so-called "F-theorem" [24]. The analyses in [23] has been extended to J\f = 2 
theories in [55] [57] [55] [IH] [SI] • See also [33] for another approximation scheme in a similar spirit. 

In the 't Hooft limit, one may deal with the planar resolvent, as in [18) . The planar resolvent is a 
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standard machinery in the context of the matrix model. It contains all the information on the planar limit 
of the matrix model. The localization formula in |16] is given in terms of a finite-dimensional integral. 
Since it looks like the partition function of a matrix model, the techniques developed in the matrix model 
context may be applicable. In jl8j . the planar resolvent was obtained for arbitrary value of the 't Hooft 
coupling, and therefore, it is possible to obtain, for example, an interpolating function for the Wilson 
loop between the weak coupling region and the strong coupling region. Systematic expansions around 
the weak coupling limit and the strong coupling limit are also available. However, it is in general difficult 
to obtain the explicit formula of the planar resolvent for generic Chern-Simons-matter theories. In the 
case of ABJM theory, a relation to a topological string theory was efficiently used, but it seems to be an 
accidental relation since the planar resolvent of ABJM theory coupled to fundamental matters does not 
seem to allow such a relation [33] . 

On the other hand, the technique in |23j is rather simple, and it can be applied to various Chern- 
Simons-matter theories. One may derive the large results directly from the matrix model action, 
without relying on the use of the planar resolvent. Up to some ansatz on the eigenvalue distributions, 
which is supported by numerical calculations, the technique provides us with the configuration of the 
eigenvalues as well as the density distribution of the eigenvalues. Therefore, following [23]: one may 
obtain all the information on the large A^ limit of the theory. However, this technique allows one to 
obtain only the results in the large A^ limit, and the systematic \/N expansion is not yet available. 

In this paper, we present the third method to investigate matrix models related to Chern-Simons- 
matter theories. It is a first-principle calculation based on the planar resolvent, but it does not require the 
expression for the resolvent as explicit as in [18]. It was shown in [34] that the planar resolvent of some 
Chern-Simons-matter theories can be written as contour integrals explicitly up to a few parameters. The 
equations which determine the parameters were also derived j34j . We show that the integral representation 
of the planar resolvent is already useful enough to derive the results in the large 't Hooft coupling limit. 
Moreover, it turns out that the derivation of those results is rather elementary; no use of the machinery 
of topological string theory is necessary. The systematic expansion in terms of an inverse power of the 
't Hooft coupling would be, though technically involved, quite straightforward. What is necessary for the 
calculation of the sub-leading terms consists of an asymptotic expansion of a given integral depending on 
a large parameter, and therefore any extra information is needed for the calculation. We demonstrate that 
such a calculation is possible for a simple Chern-Simons-matter theory. Since our method is based on the 
planar resolvent, the interpolation between the weak coupling and the strong coupling is straightforward, 
at least qualitatively. The interpolation is realized by a continuous change of the parameters, which are 
functions of the 't Hooft coupling, in the planar resolvent. 

According to our method, we reproduce known results on the large 't Hooft coupling behavior of the 
eigenvalue distributions and the Wilson loops. The theories investigated in this paper are pure Chern- 
Simons theory, U(A)fe Chern-Simons theory coupled to two adjoint matters, ABJM theory and GT theory 
all of which have at least A/" = 3 supersymmetry so that the formula in [16] is available. Here GT-theory 
refers to one of the Chern-Simons-matter theories constructed in |35| which has A/" = 3 supersymmetry 
and has the same field content as ABJM theory but the two Chern-Simons levels can be arbitrary. The 
planar limit of this theory was discussed previously in |34j . 
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This paper is organized as follows. Section [5] contains a brief summary of the theories we discuss in 
this paper. Pure Chern-Simons theory is discussed in section [3] as a warm-up. In sectional we discuss 
\J{N)k Chern-Simons theory coupled to two adjoint matters. This would be the simplest example which 
exhibits a seemingly typical behavior of Chern-Simons-matter theories with a non- vanishing sum of the 
levels. The famous results on ABJM theory are reproduced in section [5] with the method developed in 
the previous sections. In section [BJ we investigate GT theory in a manner parallel to the ABJM theory. 
It turns out that the structure of the planar resolvents of ABJM theory and GT theory are quite similar, 
but the behavior in the large 't Hooft coupling limit is qualitatively different. In section[7l the sub-leading 
order calculation is presented for the theory discussed in section S) Section [5] is devoted to discussion. 
Details of the calculations for those theories are summarized in Appendix Appendix [B] discusses the 
supersymmetry breaking in the context of the matrix model; it looks difficult to know the presence of 
the supersymmetry breaking in Chern-Simons-matter theories from the calculations in the corresponding 
matrix models. 



2 J\f = 3 Chern-Simons-matter theories and matrix models 



The family of A/" = 3 Chern-Simons-matter theories is a good arena for the study of various aspects of 
Chern-Simons-matter theories, especially in the regime where the perturbative calculation is not available. 
The Af ~ 3 supersymmetry allows one to include any number of matter fields in any representations of the 
gauge group, as long as they form A/" = 4 hypermultiplets. See e.g. [Ml- This means that there are a large 
number of theories of various kinds. On the other hand, the A/" = 3 supersymmetry is powerful enough 
to constrain quantum corrections, so that the analysis of such theories are far easier than the analysis of 
M < 2 theories. For example, the R-chargcs of the fundamental fields are not rcnormalized, and therefore 
it is rather straightforward to check whether a theory is conformal even quantum mechanically |36j . 

Following [16], we shall study J\f ~ 3 Chern-Simons-matter theories by localization technique. For a 
Chern-Simons-matter theory defined on with gauge group Y[i U{Ni)ki, the partition function is given 

by 

(2.1) 



where 



5'vcctor = - log 

I il <jl 

'5'matter — ^ ^ 5'mattcr,i? ; 
R 



sinh 



2 Ulji 



(2.2) 

(2.3) 
(2.4) 



which is applicable to any A/" = 3 theories. The explicit form of 5'inatter,_R depends on the representation 
R of the corresponding matter field. For the fundamental and the adjoint matters for U(A';) and the 
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bi- fundamental matter for \]{Ni)xlJ{Nir), it is given as 

'S'mattGrJund = ^ log 



cosh — ^ 
2 



^matter, adj 



^matter.bi— fund 



il<jl 

E 



cosh 



loa 



cosh 



(2.5) 
(2.6) 
(2.7) 



These formulae were extended to A/" = 2 theories in [53] [5S] . 

One remarkable property of the above localization formula is that the partition function (|2.ip looks 
like that of a multi-matrix model in which the angular variables are integrated out. In fact, in the case of 
ABJM theory, the partition function is directly related to a known matrix model j37j which was solved 
exactly [37] [38]. Even for other theories, the simplicity of the integrand in (j2.1|) may allow one to deal 
with them using techniques developed in the context of matrix models. 

Some quantities can be evaluated with the help of this localization formula. One such quantity is the 
free energy which was evaluated for various theories recently [22][23][26][27][28][29][31][30][32]. Another 
quantity is the expectation value of BPS Wilson loops [35]. For each unitary factor U(iV;)fe, , there is a 
Wilson loop operator Wi which preserves 1/3 of the A/" = 3 supersymmetry. The expectation value is 
given as 




(Wi) = 

where the right-hand side is defined as 

(C)mm := J Y[duisi exp[-Strcc- 



(2.8) 



(2.9) 



It is possible to discuss the 't Hooft limit, in which the ranks and the levels of every unitary factors 
U(A^/)fci are sent to infinity, while the ratios 

N, 

Xi 7^ (2.10) 

are kept finite. In this limit, the saddle-point approximation for the M-integrals becomes exact. Recall that 
such a limit can be taken for theories which does not have matter fields in higher representations. In fact, 
the allowed representations are the (anit-)fundamental, bi-fundamental, adjoint and (anti-)symmetric. If 
necessary, it is possible to assign a kind of quiver diagram, including unoriented lines corresponding to 
(anti-)symmctric matters, to each theory which admits the 't Hooft limit. Since the superpotential is 
determined completely by the supersymmetry, the quiver diagram uniquely specifies the theory. Since the 
't Hooft expansion is naturally related to the genus expansion of a dual string theory, it may be natural 
to restrict ourselves to such a sub-family of the A/" = 3 theories as long as our interest is on the relation to 
Type IIA string theory via AdS/CFT correspondence. When one would like to take the M-theory limit, 
according to [23], it might be better to consider a broader sub- family of the Af = 3 theories. 
In a suitable limit, (|2.8p can be written as 

(Wi) = / dxpiix)e^, (2.11) 
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where pi is a smooth function of compact support. Let be the supremum of the support of pi. Since 
pi is defined such that it is positive and its integral is 1, an inequahty 



{Wi) < e 



(2.12) 



holds. This implies that, if is small, then (Wi) cannot be large. Therefore, if one would like to have a 
large value of (Wi), which may be expected for a large 't Hooft coupling from AdS/CFT correspondence, 
then a;* must be large. In this case, the expectation value of Wi can be estimated as follows 



It will turn out later that the infimum of the support of pi may be negative and large when is large, 
implying that the eigenvalues are widely distributed. The arguments above is valid even when ui^t^ are 
complex, and therefore (defined in a suitable manner) is also complex, as long as Re(x*) is large. 

Although the integrand in (j2.ip are written in terms of the elementary functions, it is still difficult 
to solve the matrix model exactly, with some exceptions including ABJM theory and pure Chern-Simons 
theory. Pcrturbativc results in terms of the small 't Hooft couplings can be derived rather easily. In 
fact, one does not need to know the exact planar resolvent for this pcrturbativc calculations. This was 
demonstrated in [T7] for the evaluation of the Wilson loop in ABJM theory. 

It also turned out that the analysis in the large 't Hooft coupling regime is indeed tractable. As 
observed in [17] . a drastic simplification occurs when the width of the eigenvalue distribution is huge. 
For example, if |a;| is large, then one can use 



to simplify the saddle-point equations. According to the argument above, this approximation is valid for 
a large 't Hooft coupling. This kind of approximation is available since the complicated one-loop terms 
are written in terms of the exponential functions. A similar approximation is valid for the M-theory limit 
It was efficiently used, at the level of the integrand, to obtain the eigenvalue distribution functions in 
addition to the configurations of the condensed eigenvalues. In fact, those results contain all information 
on operators which preserve the supercharge used in the localization. As a result, the analysis of [23) 
nicely derived various results which were expected from AdS/CFT correspondence. 

It turned out [34] that there are some Chern-Simons-matter theories, in addition to ABJM theory 
and pure Chern-Simons theory, whose planar resolvent can be obtained in terms of contour integrals. 
Although the planar resolvent is not as explicit as in the case for ABJM theory, it is determined up 
to a few parameters, and the equations which determine those parameters were also obtained. In the 
following, we show that a simplification similar to the one mentioned above also occurs in the calculations 
using the integral representation of the planar resolvent, so that it is possible to reproduce some exact 
results obtained so far. One typical example of the approximation used in the following is 



(2.13) 



C01 



>th(x) = sgn(x) + 0(6-1^1) 



(2.14) 



U 



e-")(e' 



e") 




±i, {Ur <C —Or) 

^ {^ar <. Ur <. Ur) 



(2.15) 
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where Ur = Re(u) and ar — Kc{a). Here > is assumed. 



3 Pure Chern-Simons theory 

A/" = 3 pure Chcrn-Simons theory is the simplest one among the Chern-Simons-matter theories of interest. 
The matrix model [39] related to this theory has also the simplest structure. In fact, it can be solved 
easily in the large N limit [37] [35] • It is reasonable to expect that analyzing this simple matrix model will 
be very helpful to understand some structures which also appear in more complicated matrix models. 

We are interested in the large N limit. All the information in this limit can be obtained by solving 
the saddle-point equations. In this section, we consider the following equations [33] 

^u. = E-th^. (3.1) 

The indices i^j run from 1 to N. The variables Ui are assumed to be real when k is real. If k is replaced 
with —ik, then the equations (|3.ip are the saddle-point equations derived from the localization formula 
(|2.ip for JV = 3 lJ{N)k pure Chern-Simons theory. Due to this analytic continuation, it will turn out 
that Ui become imaginary. The expectation value (W) of the BPS Wilson loop [36] is given as 



W - (^E^"^ ■ (3-2) 



It is convenient to introduce new variables Zi := e"' . In terms of them, the equations p.ip can be 
written as 

log... ^ t + 2i.lE-^ (3.3) 

in the large TV limit, where 

2ttN , , 

t := (3.4) 

is the 't Hooft coupling. Following [37] [38], the planar resolvent is defined as 

v{z) ■= t dx p{x) , (3-5) 



where p{x) is formally defined as 

1 ^ 

P{x) ■■= —^5{x-z^). (3.6) 



As usual, p{x) is assumed to become a smooth function in the large N limit. 

The equation (|3.3p can be interpreted as the requirement of the balance between an external force 
(left-hand side with the opposite sign) and a repulsive force between eigenvalues (right-hand side including 
the constant term t) . Note that the repulsive force is a long-range force which is non- vanishing even for 
infinitely separated pairs of eigenvalues. As long as t is small, the eigenvalue distribution is determined by 
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the external force. It is easy to see that the eigenvalues are distributed around z = 1 at which the external 
force vanishes. Therefore, it is natural to assume that the support of p{x) is [a, b] where < a < 1 < 6 
for a non-zero t. The equation p.3p then determines the planar resolvent to be 



2tt z-x ^|(a;-a)(a;-5)| 

for z G C\[a, b]. The explicit expression after performing the contour integration is found in |37)|38). 
The definition p.Sp of the planar resolvent v{z) implies 



-t, {z = 0) 
O(z-i). (z-^oo) 



= { ^, \ r ~ \ (3.8) 



These two conditions are equivalent if a6 = 1 which is assumed in the following. Note that this is expected 
from the invariance of the equations p.ip under the simultaneous flip of the sign of Ui. The condition at 
z = implies 

dx log X 

t = ^ — (3.9) 



x^\{x-a){x-b)\ 
This equation determines a in terms of t. 

To obtain results for the pure Chern-Simons theory for which t is purely imaginary, one may start 
with this expression for a real a, and then analytically continue a, while keeping a& = 1, such that t is 
purely imaginary. Then, the resolvent p.7p with a so determined will have all the information on pure 
Chern-Simons theory in the large N limit. 

In fact, it is not difficult to perform the integration in (|3.9p . As a result, t and a := log 6 are related 

as 

e^" +e"l" ^ 2e^*. (3.10) 

This relation defines a holomorphic map from a G C to t e C. This map is a composition of the 
exponential map and the Zhukovski transformation which is well-known in fluid dynamics. If a 2> 1, 

then t behaves as 

t{a) = a-21og2 + 0(e-"). (3.11) 

By rotating the phase, a e'^a, one can make t complex. The above expression for t is still valid as 
long as \9\ < ^. However, it is not possible to deduce from p. lip the relation between t and a when t is 
purely imaginary. There is an alternative way to obtain a purely imaginary t. One may notice from the 
exact relation p.lOp that 

t{a + 27ri) = t{a) + 27ri (3.12) 

holds. Therefore, a desired value of purely imaginary t may be obtained by choosing a suitable value of 
a > which is of order one and then shifting a in the imaginary direction. This means that it is not 
enough to just specify the value of a appearing in the resolvent (|3.7p to obtain a large imaginary part of 
t. 



In the following, we show that the above results can be obtained without invoking the explicit relation 
()3.10p . For most of Chern-Simons- matter theories, it seems to be difficult to obtain an explicit expression 
for the planar resolvent, and therefore, an explicit expression for the 't Hooft coupling like (|3.10p . To 
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Figure 1: The integration contour (red line) and Figure 2: The integration contours and the 
the branch cuts (blue lines) before the shift of a. branch cuts after the shift of by 27ri. 



analyze such theories, one has to develop a technique which does not rely on any explicit expressions. 
Fortunately, at least integral representations like (|3.7p(|3.9[) can be obtained for some Chern-Simons- 
matter theories. The following calculations on the pure Chern-Simons matrix model will provide us with 
some experiences on handling those integral representations which will be applied to more complicated 
theories in the later sections. 

First, we determine the relation between t and a, assuming that the length of the branch cut is long. 
Let us focus on the case of real t so that a is also real. In terms of a new variable u := logcc, 

/'^'^ du ue^~ ,„ _ „, 
=. (3.13) 
-a ^(l-e-("+"))(l-e-("-")) 

One possible approximation valid for large a is to replace the denominator 

y^(l - e-("+«))(l - e-("-")) (3.14) 

with 1 for most of the range of li, as mentioned at the end of section[2l In fact, although this approximation 
is valid in some cases discussed later, it is not the case here. This approximation is valid if most of the 
values of u contribute equally, while here the dominant contribution is localized at u = a. A crude 
estimate of the integral p.l3p . taking this fact into account, provides 

t = a + 0{l). (3.15) 

The details of this estimate is shown in Appendix lA.ll 

The property (|3.12|) can be also derived from p.9|) . To show this, it is convenient to rewrite it as 
follows, 

dx log X , . . , 



where the contour C runs from a to 6 below the branch cut, as depicted in Figure [TJ The shift of a in the 
imaginary direction rotates the phase of a. After the shift by 27ri, the contour and the cuts are deformed 
as depicted in Figure [21 The integration along Ci gives 

dx logx f dx logx f dx 2Tri 



Ci x^y{x - a){x - b) Jc x^/{x^^^l){x^l^ Jc xy^{x - a){x - b) 

= t{a) - 2m. (3.17) 

The integration along C3 gives 

dx logx f dx logx f dx 27ri 



C3 x^(x - a)(x - b) Jci Xa/(x - a)(x - b) JCa ""^ x^(x - a)(x - b) 

dx log X 



C2 xa/(x - a)(x - 6) 



47ri. (3.18) 



These calculations show that p.l2p holds. 

To see the analytic dependence of <: on a in general, another expression 

t = ^f ^ ^"g" (3.19) 

Jc" 2m x^(x - a)(x - b) 

where C encircles the cut [a, 6], is useful. As long as a 7^ 0, there is a suitable choice of C such that 
the integrand is bounded in a neighborhood of C" and a. Therefore, the derivative with respect to a 
commutes with the integral, implying that t is analytic in a for a suitable region in C. Note that some 
singularity may appear in the limit a — > in which the contour C" will be pinched by two branch points. 
This corresponds to the divergence of t in the limit a — > 00, as seen above. This may be regarded as 
another indication that a large 't Hooft coupling is obtained only when the branch cut becomes infinitely 
long. 

The same argument can be applied to more general Chern-Simons-matter theories, as long as the 
't Hooft coupling is given in terms of contour integrals. In later sections, we understand that the 't Hooft 
coupling is an analytic function of parameters. 

The relation p.l5|) implies that, as long as Re(Q;) is large, the BPS Wilson loop behaves as 

{W) - e*. (3.20) 

As mentioned above, in the Chern-Simons matrix model, there exists a long-range repulsive force between 
eigenvalues which makes the eigenvalue distribution wider than that in the matrix models with 1 /x-type 
interaction. A matrix model of the latter kind is the Gaussian matrix model whose saddle-point equations 
are 

As is well-known in the context of AdS/CFT correspondence, the "Wilson loop" defined similarly to the 
right-hand side of p.2p behaves as 

{W) - e^. (3.22) 

This is a result showing that the presence of a long-range force may change the behavior of the Wilson 
loop for a large 't Hooft coupling in a drastic manner. 
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On the other hand, for a purely imaginary t, the Wilson loop is not larger than any function of the 
form exp[ci'''] for any c, 7 > 0. This is consistent with the exact result on the Wilson loop [40] in the pure 
Chcrn-Simons theorjifl. 



4 U(A^)fc Chern-Simons theory coupled to two adjoints 



The next equations we consider are 

k 

2^ 



— = ^coth^^-^tanh^^. (4.1) 



The indices i,j run from 1 to N. If k is replaced with —ik, then this equations arc the same as the 
saddle-point equations oiM ~ i \}{N)k Chern-Simons theory coupled to two adjoint matters (i.e. matters 
including two A/" = 4 adjoint hypermultiplets) . The number of matters is chosen such that the long-range 
repulsive force is absent. These equations were analyzed in [17] using a technique similar to [23] . 
In terms of zt = e"' , the equations (|4.ip can be written as 

logz, = 2t--y -2t--y , (4.2) 

Nj^^z,-z, N^z.-{-z,) ^ ' 

in the large N limit. This is equivalent to 

-M-(-,„ . ..lx:^-..lg^_izaL_. ,,3, 

These two sets of equations are quite similar to the saddle-point equations of ABJM matrix model [T8] . 
The main difference is that, although the planar resolvent for (|4.ip has two cuts, the positions of them 
are correlated to each other. 

The planar resolvent v{z) of this system is defined as 

viz) := t dxp{x)— t / dx p{-x)—^. (4.4) 

J Z — X J Z — X 

The saddle-point equations determine v{z) to be 



= 2. / ' (4.5) 



27r z2 - a;2 ^pT3^2y(^J— pj] 
The planar resolvent must satisfy the following conditions 



0, {z = 0) 
O(z-i). (z^oo) 



«W = { ' \ (4-6) 



^ It should be noted that the results of this matrix model for t = 2-Ki\ with A > 1 might have nothing to do with = 3 
pure Chern-Simons theory since the supersymmetry of the pure Chcrn-Simons theory is spontaneously broken when A > 1 
I41l|42ll43l . and therefore, the relation to the matrix model 113.111 is not obvious. 
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The first condition is trivially satisfied. The second one implies 



^ dx log X „ , , „, 

, 7:, = 0- (4-7) 

/a 2^^|(a;2-a2)(x2-62)| 

This is satisfied if afe = 1. In the following, we employ this choice. 
To find the relation between t and a = log 6, we have to use 

f dz v{z) , , 

t - f T- — ' (4.8) 
Jc 27ri z 

where the contour C encircles the cut [a, b], since here v{0) is not related to t. This can be written as 



i = 4 / — / — loga; , = P— 

' 2ttJ, 2tt ^ ^|(x2-a2)(x2-52)| ^2 - y2 



. dv du J {I - e-2"(i+«) 1 - e-2"(i-") 

= 2a^ — —u ^; ^ PcothQM-t; , 4.9) 

7_i 2^i_i 2t: ^(l-e-2"(i+«))(l -e-2"(i-")) ^ ^ 

where au = logx and = logy. The leading behavior of t for a large a is 

t ^ 2a / — / — u sgn(w — v) 



2a^ 

This implies that the Wilson loop behaves as 



2tt 2tt 

(4.10) 



/S7r2 \ 3 

log(W^) ^ i^—tj . (4.11) 

This coincides with the result obtained in [T7]. The behavior of (W) for a large 't Hooft coupling is 
different from the one observed in A/" = 4 super Yang-Mills theory and ABJM theory, but still it is 
exponentially increasing with the 't Hooft coupling. A crude estimate shows that sub-leading terms are 
smaller than a'^'^'^ for any e > 0. Probably, the sub-leading term is of order 0{a^). 

The analytic continuation to a purely imaginary t is now straightforward. It can be done by rotating 
the phase of a by ^. During the rotation, Re(Q;) is kept large, and therefore, the above estimate of t is 
still valid. It is also possible to perform the analytic continuation as was done for the pure Chern-Simons 
matrix model. The 't Hooft coupling t can be written as 



dy f dx log a; ^(y2 - a2)(y2 _ 52) 



/ci Jc, 2Tri - x'^ ^^^^ ^ a^){x^ ~ P)' 

where the contours Ci, C2 are depicted in Figure [31 A shift of a by deforms the cuts and contours as 
in Figure 131 For large a, the shifted quantity t{a + ^i) is estimated to be 

t(a + ^i) -t{a) - -a^. (4.13) 

V 2 / TT 

Details of the estimate are shown in Appendix IA.2I This is consistent with (14.101) . 

It is interesting to estimate the free energy of the theory. The free energy is given as 

F - /fc^u,2 + (one-loop) \ . (4.14) 
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Figure 3: The integration contour (red and 
green lines) and the branch cuts (blue lines) 
before the shift of a. 



Figure 4: The integration contours and the 
branch cuts after the shift of a by 



Assuming that the Gaussian part and the one-loop part provide contributions of the same order, it scales 
as 

F - k^Ni (4.15) 

up to an overall constant independent of N and k. This is the same scaling found in |28] in the M- 
theory limit. Since the M-theory limit in [28] is not the same as the large TV limit here, it might not be 
necessary for those results to match. However, it was shown [12] [IS] that, in the case of ABJM theory, 
the functional forms of the free energy in the M-theory limit and the 't Hooft limit coincide with each 
other. It may be expected that the behavior of the eigenvalue distribution found above would have some 
general nature for theories whose Chern-Simons levels does not sum to zero. Later, it will be shown that 
the same behavior do appear in GT theory. 

It may be interesting to see the effect of a long-range force by comparing the above results with those 

of 

A., = ^coth^+^tanh^. (4.16) 

In this system, a long-range force exists. In fact, these equations can be written as 

—Ui = ^ coth(Mi — Uj). (4-17) 

This is equivalent to the saddle-point equations (j3.ip of pure Chern-Simons theory. As was shown in 
section [21 the behavior of the eigenvalue distribution of pure Chern-Simons theory is quite different from 
that of the Chern-Simons-matter theory discussed in this section. In particular, the solution of the 
equations (|4.16|) does not have a long branch cut for a purely imaginary t, even when its absolute value 
is large. 
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5 ABJM theory 



The third example is ABJM theory. The goal of this section is to derive the famous result [TB] 

{W) - e"^, (A ^ +00) (5.1) 
where X:— from the integral representation of the planar resolvent |34] . 

Recall the saddle-point equations of ABJ theory [44] . 

— = 2^coth 2^tanh — , (5.2) 

j^i a 

~^—Va — coth /"^ _ tanh "° ""^ (5 3') 

The indices i,j run from 1 to iVi and a, 6 from 1 to A^2- Notice that there are two sets of eigenvalues, 
{ui] and {va}, according to the two gauge group factors U(A^i)fcxU(A^2)-fe- Introducing new variables 
Zi := e"' and Wa '■= e^°, the above equations can be written as 

1 \ ^ Zj 1 \ ^ ( — w„) 

logz, = ti - ^2 + 2ti ■ — ^ -2t2-—J2 \ V 

\og{-{-Wa)) = t,-t2 + 2t,-±Y.r^, 2t2.-l^ (5.5) 

A^l ^ (-tWa) - -Zi ^2 (-t"a) - ("Wb) 

where 

h ■■= h (5.6) 

In the following, ti and ^2 are regarded as real variables, and they will be analytically continued to the 
above values later. The structure of the saddle-point equations suggests that there would be two cuts, 
one is around z = -f 1 where Zi condense, and the other is around z = —1 where —Wa condense. Define 
the planar resolvent 

dx p{x) t2 I dx p{x) . (5-7) 

Z — X J Z — X 

The distribution of Zi is described by p{x), and that of —Wa is described by p{x). We assume that the 
supports of p{x) and p{x) are [a,b] and [c, d], respectively, where & < < c. The saddle-point equations 
determine v{z) to be 



log(-e*2-*i2:) ^(z - a){z - b){z - c){z - d) 



v{z) 



2tt z-x y^\{x - a){x - b){x - c){x - d)\ 
The conditions imposed on v{z) are 



a 27r z-x y/\{x - a){x - b){x - c)ix - d)\ 
dxlog{e*^'^^x) y/{z - a){z - b){z ~ c){z - d) 



(5.8) 



t2-ii, (2 = 0), 
O(z-i). (z^oo) 



= \ -1. , X (5-9) 
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This amounts to three conditions on the parameters a, b, c, d. These are reduced to one condition, assum- 
ing 

ab =1, cd = 1. (5.10) 

The remaining equation 

t2-ti _ _ dx log(-a:) 

2 ^ 27r xy/\{x - a){x - b){x - c){x - d)\ 

+ f ^ , ^"g" (5.11) 

Jc 2-K x^\{x - a){x - b){x - c)ix - d)\ 

provides a relation among the undetermined parameters and the 't Hooft couphngs. One more condition 

is necessary to completely determine the parameters. One may use one of the following relations 

f dz viz) f dz v(z) 

tl = i TT- — . h = - f ^ — , (5-12) 
Jc^ 2tti z Jc_ 2m z 

where C+ encircles the cut [c, d] and C_ encircles [a, 6]. It was shown in [34] that the relation (jS.lip can 
be written explicitly as follows, 

{c + d)-{a + b) = 4e*i-*^ (5.13) 



First, let us focus on the case ti = t2, that is, the ranks of two gauge groups are equal, Ni = N2 = N. 
Suppose that both \a\ and d are large for which the corresponding 't Hooft couplings are large. Then, 
the relation (|5.13[) implies 

a + 5 = A-2, c + d = A + 2, (5.14) 

where |A| is large. This relation can be satisfied only after an analytic continuation, since originally 
a < < d was assumed. The resulting branch cuts should be almost parallel to each other, but the 
direction of the cuts in the complex plane is not determined by (|5.14|) . 

Here, recalling the exact solution found in [18] will be helpful to find a more detailed information of 
the branch cuts. It was shown that 

a = 2 + iK, (3 = 2-iK (5.15) 

holds, where a = c + d and (3 = —a ~ b in our notation. Here k is a real function of the 't Hooft coupling 
which is large and positive when the 't Hooft coupling is large. This indicates that the lengths of the 
cuts are huge as expected above, and the cuts are almost parallel to the imaginary axis. 

With this input, we show that the integral representation of the planar resolvent ()5.8|) can reproduce 
the relation between the 't Hooft coupling and the length of the branch cuts. We assume 

a = -ie'> - 2 + 0{e-'>), d = -ie'> + 2 + Oie"*), (5.16) 

where 7 > is assumed to be large. Of course, this is compatible with (|5.14p . It can be shown that the 
flip of signs in the imaginary part of a, d results in the flip of the sign of the 't Hooft coupling. 
The 't Hooft coupling ti can be written as 



U = 2 



ci 



dy 
2111 



dx 


log a; 




a){y- 


b){y- 


c){y- 


-d) 


2ni 


y{y~x) 




a){x — 


b)ix~ 


c)(x 


-d) 


dx 


1 




-o){y- 


-b){y- 


-c){y 


-d) 



+2m I ^ I ^ , ^ , V " ^'"^ ^' . (5.17) 

Jo, 27ri Jci 2iTi y{y - x) y/{x - a){x - b){x - c){x - d) 
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Figure 5: The integration contours (red and Figure 6: The deformed integration contours 

green lines) and the branch cuts (blue lines) convenient for estimation, 

before the shift of a. 



The contours are depicted in Figure [S] To estimate the integrals, it is convenient to deform the contour 
as in Figure |6l It turns out that the first term in the right-hand side is of order 0(7). The dominant 
contribution coming from the second term is 



Jet 27ri Jc- 2Tn y{y ~ x) ^{x - a){x - b){x - c){x - d) 

The other terms are negligible compared to this term. By estimating this integral, one obtains the 
following asymptotic behavior 

h ~ -7'- (5.19) 

IT 

See Appendix IA.3I for the details. Recall that ti = 27riA and A = Then 7 can be solved in terms of A 
as 

7 ~ ttV2X. (5.20) 



This then implies the desired behavior (|5.ip of the Wilson loop for large A including the exact coefficient 
in the exponent. 

In retrospect, it would be possible to find the right configuration of the branch cuts without relying on 
the exact solution in [TS]. One may start with the ansatz (|5.16p which can be deduced from (|5.13p . One 
may try with a real 7, and obtain the result (|5.19|) . Recall that (|5.13|) allows a complex 7. Since t should 
be purely imaginary, one cannot rotate the phase of 7. It is still possible to shift 7 in the imaginary 
direction by an 0(1) amount. This does not change (|5.19p . but docs change the 0(7) term by a real 
amount. Therefore, whether such a shift is necessary or not can be determined by estimating the 0(7) 
terms. This then determines the direction of the cuts. Anyway, if one is only interested in the leading 
order behavior, the direction of the cuts turned out to be irrelevant. 
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The possibility of an interpolation between the weak coupling region and the strong coupling region 
for. say, the Wilson loop is obvious. As was shown in [T7] that the small 't Hooft coupling limit corresponds 
to a limit in which the branch cuts of the planar resolvent shrink to a point. It is easy to imagine a 
continuous deformation of the branch cuts from the point-like one to the one depicted in Figure [Sj while 
keeping the conditions ab = cd = 1. This indicates that the weak coupling results and the strong coupling 
results can be connected continuously. 

Next, let us consider the case ti ^ t2 and both ti and t2 are purely imaginary. Recall the relation 

{c + d)-{a + b) = 4e*i-*^ (5.21) 

The right-hand side is of order one, and therefore, the positions of the cuts arc the same as those of 
ABJM theory at the leading order of the large 't Hooft couplings. A natural guess for the sub-leading 
terms would be 

a ^e''-2e*l-*^ d ie'' -f 2e*i-*^ (5.22) 

It can be shown that the resulting relation between 7 and ti are the same at the leading order of 7. As 
a result, the Wilson loops in ABJ theory behaves as 

{W) ^ exp [ttVAi + A2] , A1.2 := (5.23) 

This was obtained in |18| . 

It was claimed in [44] that ABJ theory is well-defined only for \Ni — < k. This implies 

|ti-i2|<27r (5.24) 

is allowed to consider. From the matrix model point of view, however, it seems to be difficult to find any 
sign of the ill-definedness for the corresponding parameter region. Indeed, an analytic continuation by 
which (|5.24[) is violated does not show any singular behavior in the planar resolvent, nor in the Wilson 
loop. It is similar in pure Chern-Simons theory where the matrix model is well-defined for any value of 
the 't Hooft coupling. In addition, the Wilson loop docs not show any special behavior at A = 1. In 
the appendix [Bl the analysis on A/" = 3 \i{N)k Chern-Simons theory coupled to fundamental matters is 
summarized. In this theory, the supersymmetry is broken for a choice of parameters which is expected 
from the argument based on the brane construction. However, there is also no special behavior of the 
Wilson loop at a value of the 't Hooft coupling at which the supersymmetry breaking is expected to 
occur. 



6 GT theory 

One of the straightforward generalization of ABJM theory, at least from the point of view of saddle-point 
equations, is GT theory [35]. This is proposed to be dual to a massive Type HA string theory [35] [45]. 
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The localization procedure of |16| can be applied to GT theory. The resulting saddle-point equations are 

- — Ui = > coth — > tanh , (6.1) 

■^^Vn = coth -2 ^ — y^ tanh -2 -. (6.2) 

27ri ^ 2 ^ 2 ^ ^ 

When fci + A:2 = 0, GT theory is reduced to ABJM theory. The planar resolvent defined as in (|5.7p was 
obtained in [SI] in terms of the following contour integrals 

/■'' dx log(-e"^^*^~*^'.T) J{z - a){z - b){z - c){z - d) 

viz) = K2 

A 2^ z-^ ^\ix - a){x ^ b){x ^ c){x - d)\ 

f' dx log(e^'*--*^)x) ^{z - a){z - b){z - c){z - d) 

+Ki / , =, (6.3) 

7c 2^ z-x ^\{x - a){x ^ b)ix - c)ix - d)\ 

Ki,2 := — , ^ . (6.4) 

All the parameters fci^2,-/Vi,2 are assumed to be proportional to a common number k which is sent to 
infinity. As in the case of ABJM theory, assuming ah ^ 1 and cd = 1, the number of the parameters in 
the planar resolvent can be reduced to two, and the remaining parameters are then related to the 't Hooft 
couplings ti^2 as (|5.12p . Instead, one of these relations can be replaced with 

t2 - h dx log(-a:;) 

— 7i — = ^2 



where 



a 27r x^y\{x - a){x ~ b){x -~ c){x - d)\ 
dx losx 



+Ki / , (6.5) 

Jc ^TT xy'Kx - a){x - b){x - c){x - d)\ 

Let us consider the case where ki + K2 7^ and ti — ^2- In this case, an explicit expression for ()6.5|) like 
()5.13p is not known. Therefore, it is necessary to extract some information on the branch cuts directly 
from this integral representation. 

We assume that ki and K2 are of order one. Let a = — e" and d ~ , and both a and (3 are large. It 
turns out that an assumption /? 3> a is not compatible with (|6.5p since in this case the right-hand side 



is estimated to be ^. It is of course similar for the other assumption a 3> /?. Therefore, it should be 
assumed that a and (3 are of the same order. This assumption then implies 

K2e"-Kie'^ - 0. (6.6) 
This relation also holds for the ABJM case ki = — K2 = 1- Indeed, in this case there is one solution 

a ~ 7+^*, P ^ (6.7) 
where 7 is a large and positive number for which 

a ^ —ie^, d ^ ~ie^ . (6-8) 



This is the ansatz used in the previous section. Since ki and K2 are of order one, the solution of (|6.6 
may not be different so much from the above one. That is, the solution should be the same with (j6.8 
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up to an order one correction. Note that the solution (|6.8p is also valid for a complex 7 as long as Rc(7) 
is large. 



Now it is possible to estimate ti as a function of 7 in the large 7 limit. We assume first that 7 is real, 
and its phase will be rotated if necessary. The formula for ti is 



-2ni 



dy 



dx \ogx \J{y- 


a){y- 


b){y- 


c){y- 


d) 


2iTiy[y-x) [x - 


a){x — 


b){x- 


c){x- 


d) 


dx log X \/(y~ 


a){y- 




c){y- 


d) 


2Triy(y-x) y^{x- 


a){x — 


b){x- 


c){x - 


d) 


' dx 1 y/{y 


- a){y 




- c)iy 


-d) 



n f dy I "-"^ ^ V \ti "-/ya "i\a '-j\t) "y /n q\ 

~ Jc, Jc^ y{y - x) ^ {x ^ a){x - h){x - c){x - d) 

Since the configuration of the branch cuts are almost the same with that for ABJM theory, the estimate 
of the integrals proceeds in a similar way. Especially, one may immediately find that the third integral 
provides an 0(7^) contribution. 

The difference from the ABJM theory comes from the other two terms. One can show that the leading 
contribution to each integral is of order 0(7"^). It is larger than the third integral because of the presence 
of log a; in the integrand which is of order 0(7) in most of the integration region of x. What happened in 
ABJM theory is that this leading contributions exactly cancel between the first and the second terms, so 
that the second largest term coming from the third integral determines the asymptotic relation between <i 
and 7. Here in GT theory, such a cancellation is not complete, and the remaining 0(7'^) term dominates. 
The resulting relation is 

~ (6,10) 

See Appendix IA.4I for the details. The right-hand side is real when 7 is real. To obtain the purely 
imaginary ti, one may rotate the phase of 7 by ^. During this rotation of 7, Re(7) is kept large, and 
therefore the above relation is valid. In terms of A := 7 is given as 

7 = \i\e^\ |7l ^ (^—A ■ (6-11) 

Note that in (|6.11|) the number k cancels out, so 7 is independent of the choice of fc, as it should be. 

This result is qualitatively similar to the result obtained in section S) As discussed there, this scaling 
suggests that the free energy of GT theory scales like fc^TVi which coincides with the one found in [28] 
although the limit taken there is not the usual 't Hooft limit considered here. Interestingly, not only the 
scaling but also the numerical factor is the same with the result in [28] . In [28] , the eigenvalues lie on the 
line 

y = (6.12) 

which corresponds to the analytic continuation of 7 mentioned above. The maximum value of x is 
determined by p{x^,) = where 

27rp3 3-\/37r^p' 
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where k := ki + k2. The necklace quiver with p ~ 2 in |28] is GT theory. One obtains 

which coincides with the result obtained above. 

Another immediate consequence of (|6.1ip is that the Wilson loop of GT theory behaves as 



(6.14) 



exp 



2 \Kl + K2 



(6.15) 



Since A^i = N2 was assumed in the above calculations, this formula is valid for Wilson loops for both 
U(A^) factors. This is actually consistenjfl with the results in [57] where it was argued that the radius / 
of the dual AdS4 in massive Type IIA supergravity is 



I 



(6.16) 



where hq is proportional to the _Fo flux. Since the exponent of should be proportional to an area in 

the AdS4, the Wilson loop should scale like (|6.15|) . Note that the scaling (|6.16|) was shown to be valid if 



^3 



(6.17) 



which is indeed satisfied in our setup. 



7 Sub-leading contributions 



In this section, we present the calculation of the sub-leading contributions to t{a) for the theory considered 
in sectional To do this, it turns out that the following expression for t{a) is convenient. 



+a 



du 



2711 J _^ 27r ?/2 _ (,2u _ g-2(Q+M))(-2 _ g-2(Q-u)) ' 

where the contour C is depicted in Figure [71 

The dominant contributions come from the following parts of the integral 



du we" 



V(2/2-a2)(j^2_^2) 



dy_ 

2m J _^ 2-K y'^ - e^" ^(1 - g-2(Q+«))(-2 _ g-2(a-«)) 



27r2 



dv 



+a 



du u sgn(w — v) 



-\u—v\ 



cosh(it — v) 



v/(l + e-2(°+^))(l + e-2("=^ 
^/(l - e-2("+"))(l - e-2("-«)) ' 



(7.1) 



(7.2) 



^ We would like to thank Soo-Jong Roy for drawing our attention to this issue, and checking independently that the 
scaling I I6.15I I agrees with supergravity. 
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Figure 7: The integration contour (red and green lines) and the branch cuts (blue lines). 
This integral is divided into two parts. The first part is evaluated as follows, 



— / dv duusgn(u — v) — p= 



v/(l + e-2("+"))(l + e-2("-'')) 

- e-2("+«))(l - e-2(a-«)) 



2 r , r , v/(l + e-2("+^'))(l + e-2(a-i,)) 
— / du dvu- 



./o 



'^(1 _e-2("+«))(l -e-2("-")) 



Jo V 1 - e-2(a-") 



(7.3) 



where exponentially small terms are neglected. Now w-integration can be performed exactly. As a result, 
one obtains 



„ , dfii — 

7r2 Jq Vl - e-2(a-") 
2 



+ >/l + e-2(a-«) _ 1 + log2 - log(l + + e-2(a-«)) 



37r2 



+ 0(a), 



(7.4) 



This coincides with (|4.10p . Note that the coefficients of 0{a) terms can be determined. For example, an 
0{a) integral can be estimated as follows. 



— / du 



Vl - e-2(a-") 



(Vl + e-2("-") - 1) 



2 r , ia-t) , r- 
dt ^ ^. (Vl 



2t 



1) 



2 f°° 1 

a/ , _^^ (Vl + e-2t-l) + C>(l). 



•7r2 



vT- e- 



(7.5) 



The last equality holds since the integrals in the second line are convergent when the integration region 
[0,a] is replaced with [0,oo]. 
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The second part turns out to be of order 0{a). Indeed, 



' dv duu sgn(it — v) - 



27r2 7„„ " ^C0Sh(u - v) ^(l - e-2(a+«))(l _ e-2(a-«)) 



The contributions from the other parts of the contour C turn out to be of order 0{a). For example, 



dy du ue-" {v"^ - a'^){y^ - h"^) 



2iTi 27r ?/2 _ ^(1 - e-2("+"))(l - e-2("-")) 



27r2 y__^ '^^e2E+2»« - e2("-") ^(1 - e-2("+"))(i _ e-2("-")) 

2^ "^^io '^''e2^+2i«-e2("-") ^1 _ e-2(a-n)= ^ 



= de dt , ^ — + (7.7) 

In summary, the asymptotic expansion of t{a) for large a is 

= -^a^ + cia + 0(l), (7.8) 
where ci is a constant which can be calculated at least numerically. Inverting the relation, one obtains 

37r2\ ^ 1 1 /37r2N ^ 



«(0 = ( — ) ^'-gl — ) ^"'+0(t-^). (7.9) 



8 Discussion 



We have shown that some exact planar results on the eigenvalue distributions in matrix models related 
to Chern-Simons-matter theories can be derived from integral representations of the resolvents. All the 
results derived in this paper were already derived using a topological string theory |18j or a saddle-point 
approximation for the localization formula of the partition function 1:23^, but the method in this paper is 
new. Our method requires the determination of the planar resolvent, but it does not need to be as explicit 
as in the case of ABJM theory. In addition, the method to derived the information on the eigenvalue 
distribution is a rather elementary estimate of an ordinary single- or double integrals. The systematic 
expansion around the large 't Hooft coupling limit is possible without any conceptual difficulty. Since the 
resolvent is determined for any value of the 't Hooft couplings, up to the exact positions of the branch 
points, the existence of an interpolating function of, for example, the expectation value of the BPS Wilson 
loop {W) is evident. 
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One of the results obtained in this paper is the relation between the 't Hooft coupling t and a parameter 
a which determines the position of a branch point. In general, the relation would be of the form, 

t = c^a^ + ic2a'^ + cia + ■ ■ ■ , (8-1) 

provided that Re(Q;) > is large. The coefficients ci.2,3 arc real. In general, C3 is non-zero, and then 
t behaves as even after an analytic continuation which brings t to be purely imaginary. This seems 
to be typical when the sum of the Chern-Simons levels are non-zero. Such a behavior turns out to be 
consistent with the conjecture [35]. For a suitable choice of the theory with a suitable choice of the 
parameters, C3 = may be realized while C2 7^ 0. Then the theory exhibits some properties in the large 
't Hooft coupling limit which is expected from an AdS4 geometry in an ordinary (massless) gravity dual. 
In addition, there exists a special kind of theories for which C3 = C2 = 0. In such theories, the properties 
in the large 't Hooft coupling limit may be quite different from the other two kinds of theories mentioned 
above. It is interesting that our analysis based on the planar resolvent provides us with a rather unified 
picture on the structure of Chern-Simons- matter theories in the large 't Hooft coupling limit. 

It seems that our method can be applied to more general theories. Our method would be useful since 
it does not require one to obtain the planar resolvent in the most explicit form. In addition, it seems 
that our method requires less numbers of assumptions compared with [23j . 

It would be interesting to apply our method to a family of circular quiver Chern-Simons-matter 
theories studied recently in |23j . In their research, it was found that the eigenvalue distribution is not 
always linear but it can be piecewise linear. It seems to be difficult to understand this fact from the 
viewpoint of the planar resolvent since, usually, the eigenvalue distribution corresponds to the branch 
cuts of the resolvent, and the cuts are usually taken to be linear. Of course, the configuration of the 
branch cut is not determined a priori, and therefore, the branch cuts can be chosen to be piecewise linear. 
However, if it is the case, then it should be possible to know the reason why it should be so. It is also 
possible that the piecewise-linearity might suggest that the method based on the resolvent does not work 
for those theories. 

Assuming that our method could be applied to a large family of theories, it would be interesting to 
classify the behavior of the eigenvalue distribution in the large 't Hooft coupling limit. In this paper, 
we found three patterns, two of which have geometric interpretation via AdS/CFT correspondence. 
It would be very surprising if these three exhaust all possible behaviors. If it is the case, one might 
interpret this as an indication that the presence of a gravity dual is rather common among Af = 3 
Chern-Simons-matter theories with a sensible 't Hooft limit. If there exist other behaviors, then it would 
be interesting to investigate the implication of such behaviors, possibly related to a new example of 
AdS/CFT correspondence. 

The main task for the calculations in this paper is to obtain the asymptotic expansion of a function 
t{a) which is defined in terms of a double integrals with the integrand including logx. It would be 
interesting if there exists a mathematical framework dealing with such functions, enabling one to make 
our calculations transparent. 
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A Estimate of integrals 

A.l Pure Chern-Simons theory 

We evaluate the following integral 



This is divided into three parts /i + /2 + /a where 



J _ r du {u- a)e^("-") 

^ Jo 71" ^(1 - e-("+«))(l - e-("-")) ' 



(A.1) 



^1 / =, A.2 

Jo y'(l-e-("+"))(l-e-("-")) 



(A.3) 
(A.4) 



Noticing 

^ < ^ < (A.5) 

for u e [0,a], /i satisfies 



, ,1 sin^^e"- < /i < , - 1 sin"^e"2 , (A.6) 

Vl -e-2" V ^ / " VI -e-" V / 

which implies 

/i = a + 0(ae-^). (A.7) 

/2 can be estimated as follows. 

/a < \ y I 



\J\ - e-" 7o 7^ Vl - e"-" 
< , ' rdtte^^ 
Vl - e"" Jo Vl - e-* 

Therefore, 

/2 = 0(1). (A.9) 



24 



The remaining part I3 can be estimated as follows. 

du 



I" v^(l - e-("+"))(l - e-("-")) 



< 



ae 2 



+ - log(l + Vl - e-") 



Therefore, this part is negligible. Combining these estimates, we obtained (|3.15p . 



(A.IO) 



A. 2 Chern-Simons theory with two adjoints 

In section m an analytic continuation of 



'(») = 1/5, 



dx log X ^ — a'^){y'^ — h"^) 
Ic, 2t:i Jc, 2^y23^ yp=pyp=|fy 

is discussed, where b = e" and ab = 1. Shifting a by this becomes 

dy 



(A.ll) 



/ TT 



dx log X ^ (y^ + a?){y'^ + 6^) 



(A.12) 



where the integration contours are depicted in Figured To estimate the resulting integral, it is convenient 
to divide it into four integrals labeled by contours. 

The first integral wc estimate is the one for the contours (C[, C2). This can be written as 



= 4 



dy 



dy 



dx log X (y2 - 


- a'^){y^ - 


h62) 


2m y^ - x-' ^(x^ -\ 


- a2)(a;2 - 


f 62) 


dx log(ia;) ^ - 


- a^)(y^ - 


-62) 



C, 2TTi Jc, 2TTi y2 _ ^(2;2-a2)(a;2-62) 
dx 



t(a) + 2m I ^ 



1 



V(2/^-a2)(y2„52) 



27ri 27ri 2/2 - ^(x2 - a2)(2;2 - 62) 
The second term of the last line is estimated as follows. 



(A.13) 
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TT 
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Next, consider the contours (C(,C4). The corresponding integral can be estimated as 

4 



dy f dx \ogx \J {y^ + a?){y'^ + 



dy ^/W+~a^)WT¥) d(f)-ae-°'+"^ 



- a2)(b2 - ^2) 
27r y2 



2i 

— a 



2 



(A.15) 



It turns out that the remaining integrals are negligible compared to the above terms. Therefore, we 
obtained the estimate (|4.13p . 



A.3 ABJM theory 



It turns out that the dominant contribution to the 't Hooft coupling ti comes from the following integral 

dx 



2m 



dy 



1 



\/{v- a){y - b){y - c){y - d) 



C+ 27ri Jc- 27ri y{y - x) ^(a; - a){x - h){x - c)(x - d) 



where 



-ie^ - 2 + 0{e-^), d = -ie'' + 2 + 0{e-''), ab = cd = 1, 



(A.16) 
(A.17) 



and the integration contours are depicted in Figure [6] This can be written as follows. 

f.-ie+''+Q 



27ri 



27ri 



27ri 
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The term with the delta function provides a contribution of order 0(7). The remaining term can be 
estimated as follows. 



2m 



dy 



dx 1_ 



1 ^ {y ~ ia){y - ib){y - ic){y -id) 
--, 2m J ^+-, 2m y' y - x ^[x - ia){x - ib){x - ic){x - id) 



„ dv du ^ e" 
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A.4 GT theory 



The dominant contribution to <i in (j6.9p comes from the following two integrals 



-2-2/ ^ 



dx loj 
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c){y 
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2711 y{y ^ 
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+2-1 / , / , \r ":r ::r ":r (A.20) 



Note that they cancel each other if ki = K2. In general, it is expected that these terms provides the 
contributions of order 0(7'^); the length of the integration range in terms of logx and logy is of order 
0(7) and logcc in the integrand provides another 0(7) contribution. This 0(7^) contribution vanishes 
only for the ABJM slice, so that the deformation from ABJM theory to GT theory is not continuous in 
the large 't Hooft coupling limit. 

Each integral can be estimated similarly. Let us focus on the first one. By the similar calculation in 
the ABJM case, it can be estimated as 



I — ^ '^^ ~ - b){y - c){y - d) 



c+ 27ri Jc- 27ri y{y - x) ^{x - a){x - h){x - c){x - d) 



„ , +^ dv du ^ ue" 
2-2 / / P 

' _^ 27ri 2m e" - e" 

dv [+'' du , u 

-2 / 7: — : / 7: — : u P coth 



^ 27ri 27ri 2 

3 f^^ dv f^^ du 

37r2^ ■ 

Note that an extra sign appears in the second integral in (jA.20p due to the opposite direction of the 
integration contour for x. As a result, the 't Hooft coupling ti behaves as (|6.10p in the large 7 limit. 
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B On SUSY breaking in Chern-Simons-flavor theory 

Consider Af = 2 pure Chern-Simons theory. The expectation value of the normalized BPS Wilson loop 

is [m 

oc 4^^^^^ (B.l) 

^ ' N sin(7r/A:) ^ ^ 

up to an overall phase. 

It is easy to see that {W) is non-zero for < x ^^'^ vanishes at ^ = 1. Curiously, = 1 is 
the boundary of the parameter region in which the supersymmetry is spontaneously broken [41] |42) [43] . 
Note that this behavior of the Wilson loop is preserved in the 't Hooft limit. Indeed, 

^.^ 1 sin(7rA^/fc) ^ ^.^ 1 sin(7rA) ^ sin(7rA) 
N^oo N sin(7r/fc) N^oo N sminX/N) ttA ' ^ ' ' 



27 



which vanishes at A = 1 but is non-zero for < < 1, where A = 

This might suggest that the vanishing of the Wilson loop might be related to the breaking of super- 
symmetry. To check whether this could be the case, we consider another example, A/" = 3 Chern-Simons 
theory coupled to Nf fundamental hypermultiplets. 

First, let us consider the expected pattern of the supersymmetry breaking in M ~ 3 Chern-Simons 
theory with flavors. We assume that the supersymmetry breaking pattern of this theory would be the 
same as Af ~ 2 Yang-Mills Chern-Simons theory with the same gauge group coupled to the same number 
of fundamental hypermultiplets. (We expect that the deformation from TV = 2 theory to JV = 3 theory, 
which is a continuous change in the superpotential, does not change the vacuum structure.) Then, the 
supersymmetry breaking pattern can be deduced from the corresponding D-brane configuration. 

The brane configuration in question consists of two parallel NS5-branes, N D3-branes suspended 
between the NS5-branes, k D5-branes which are intersecting with one of the NS5-branes, and Nf D5- 
branes which are parallel to the other set of D5-branes and are intersecting with the D3-branes. Since 
the D3-branes have a finite extent in one direction, the worldvolume theory is a three-dimensional one. 
In the absence of the second set of D5-branes, the worldvolume theory on the D3-branes is N = 2 \J{N)k 
Yang-Mills Chern-Simons theory, and Nf D5-branes add massless fundamental hypermultiplets to the 
theory. The question of whether the ground state of the theory preserves the supersymmetry can be 
answered using the s-rule [IS]. Roughly speaking, this rule says that only a single D3-brane can be 
suspended between an NS5-branc and a D5-branc in the above setup. According to this rule, only up to 
k D3-branes can be suspended between the NS5-branes. This provides, in the pure Chern-Simons theory, 
the limitation < 1 for the 't Hooft coupling for which the supersymmetry is preserved. In the presence 
of Nf D5-branes, D3-branes between the NS5-branes can be broken by these D5-branes. Therefore, more 
D3-branes can be present without breaking supersymmetry. The upper bound on N is therefore 

N<k + Nf, (B.3) 

implying 

N Nf , , 

The inequality (|B.3P can be also understood in the following way. One may start with a configuration in 
which all the k + Nf D5-branes are intersecting to an NS5-brane. For this configuration, the condition 
for preserving the supersymmetry is (jB.3p . One then separates Nf D5-branes from the intersection and 
obtains the configuration discussed above. This procedure consists of the web-deformation of the D5-NS5 
system which does not break any supersymmetry. Therefore, the condition (jB.3| is the condition for 
preserving the supersymmetry for the flavored Chern-Simons theory. 

It is interesting to see whether the Wilson loop in this flavored theory vanishes when the inequality 
(jB.SP is saturated. 

The planar resolvent of A/" = 3 Chern-Simons theory with flavors was obtained in [34] . The Wilson 
loop is obtained from the resolvent. The explicit expression is as follows. 



{W) = - 



(a-a-i)2 t2 



h [A{2 + a + a-^) - 2 



-t2 + -^\/2 + a + a-^ 



(B.5) 
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where ti oc ^2 oc and a is determined by 



log- 



+ <5 



It is convenient to introduce u such that a — e^". Then we obtain 

1 



(W) 



sinh M + t2(cosh u — 1) 



sinh u + 2^2 sinh — 



where u satisfies 



2 log coshw 



cosh 1 



f -?. - U 



Note that t2 ~ Q implies coshu = e'^/^ which then implies 



ii/2 



Recalling ti = 27riA, we found that the formula (|B.7I) reproduces the known result (|B.2[) 
The vanishing of the Wilson loop requires 



,2 W *2 

cosh — = , 

2 2 ' 



or 



coshw = — 1 — t2- 
Therefore, the Wilson loop vanishes only when 



l0g(l + t2f - 



l+t2 



= t2-tl. 



(B.6) 



(B.7) 
(B.8) 

(B.9) 



(B.IO) 
(B.ll) 

(B.12) 



This has nothing to do with the supersymmetry breaking condition (jB.4[) . 

Note that the condition (|B.4|) for the supersymmetry breaking is consistent with the result in [34] . In 
[34] . we started with GT theory whose gauge group is \]{N)kiX\]{N)^k2- If = ^"2, then it is ABJM 
theory, and therefore the supersymmetry is preserved for any value of the parameters. Then, we took the 
limit fc2 oo. This limit is a weak coupling limit in terms of the second gauge group, and therefore, the 
supersymmetry should not be broken. The resulting theory is a Chern-Simons theory with gauge group 
U(A^)fcj coupled to Nf ~ 2N fundamental hypermultiplets. This theory satisfies the condition (jB.4p for 
any choice of N and ki , so the supersymmetry is preserved. This must be the case since the localization 
formula used in [34] is valid in the presence of the preserved supersymmetry. 
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